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Abstract 

The recent accurate measurements of Berg, Moldover and Zimmerli 
^D ! of the viscoelastic effect near the critical point of xenon has shown that 

^L) I the scale factor involved in the frequency scaling is about twice the 

f^ ■ scale factor obtained theoretically. We show that this discrepancy is 

^D , a consequence of using first order perturbation theory. Including two 

loop contribution goes a long way towards removing the discrepancy. 
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The shear viscosity 77 of a hquid-gas system near the critical point or a 
binary hquid mixture near the critical consolute point has a weak divergence 
characterized by a small exponent a;^. If e is the deviation of the temperature 
T from the critical temperature Tc,i-e e=(T-Tc)/T, then the correlation length 
C, of the fluctuation diverges as e ^ 0, according to .^ ~ e"" The divergence 
of the shear viscosity is expressed as 

V oc T" (1) 

If the measurements are conducted at a finite frequency , a new length scale 
enters the problem. This is due to critical slowing down , which implies 
that the time r of the fiuctuations of size k~^ (k is the wave number of 
the fiuctuation) diverges as r oc k~^ for small k. This produces a length 
scale r^/^, which is tu^^^^ where to is the frequency at which process is being 
probed. If the external frequency goes to zero, this length scale l^ oc cj^^/^ 
goes to infinity and ^ is the only controlling length in the problem. For a 
finite value of l^ ,it is possible to go also enough to the critical point ,such 
that ^ ^ luj and then as ^ exceeds l^^, the viscosity can change no longer and 
is determined by l^. Thus for ,^ ^ 00, 

r] oc e oc cj-""/" (2) 

In D-dimension ,the exponent z=D+a;^ and a finity values of C, and u ,the 
viscosity is described by the scaling law (7,8), 

where Tqk^ is the characteristic frequency associated with the decay of fiuc- 
tuations and S is a scaling function characterized by S(0)=constant and S(y) 
oc y for y>>l. The above constraints on S(y) helps us recover Equations (1) 
and (2). In view of the smallness of the exponent a;^ (~ 0.067), we can expand 

^'''' as 1+Xrjln^ and writing [S{y)]~'''^^^^'''^=l-[x^/{D + x,j)]lnS{y) + ,we 

arrive at 

Ar] = r]{^, uj) - r]{^, 0) = -r]{^, 0) [x,/{D + x,,)]ln{cj/ToK'') (4) 

The simplest scaling function that one can think of is 

lnS{n) = ln{l+an) (5) 

where fl = uj/{2rQK,^) = 2kT % 7 which is the scaled frequency with the 

frequency scale set by the Kawasaki form and "a" is a number of 0(1) which 
describes where the crossover from the "hydrodynamic" (zero frequency) to 
the " no n- hydrodynamic" (frequency limited) behaviour takes place. 

In the one-loop self consistent calculation, the frequency dependent shear 
viscosity in D=3 is given by (see Fig.l) 

1 /■ d^p Ap'^sin^Ocos'^O 

Viii^,^) = i J (2vr)3 (p2 + ^2)2 [_,^ + 2TopWp' + f^''] 

(6) 
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leading to 

Comparing with Equations (4)and (5) ,the function lnS(r2)to one loop is 
In5'i(f2), given by 



lnSi{VL) 



p^dp 1 1 . 



wliicli has the high frequency form, lnS'i(r2)~|lnr2/e'^~^'"2^ corresponding to 
"a" =0.147. Yet another way of estimating "a" is to study the low frequency 
limit of Eq.(8), from where we find 

InSiiQ) = -iQ^ + 0{Q^) (9) 
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which corresponds to "a"=||=0.589, when we compare with the low fre- 
quency Taylor expansion of Eq.(9). Thus, we see that for the crossover 
frequency scale "a", there are two possible estimates (8). One comes from 
the high frequency end, which we call ah^ and another from the low frequency 
end, which we call a^g. At one loop level 0/^^=0.147 and a;(,=0.589. The two 
estimates are quite different. The full scaling function, obtained by evaluating 
the complete integral in Eq.(8), gives the gradual change in scale from a^^ to 
Qhi- This was done in Ref.(8). 

In recent accurate viscoelastic effect measurements of Berg et al.(9,10), 
it was reported that the true frequency scale is about twice as big as the 
theoretical one. Since the one loop scaling function has a changing frequency 
scale, it is worthwhile examining this observed discrepancy more closely. Ac- 
cordingly, we analysed the data of Berg et al. slightly differently. For the 
different t= ^ " values frequencies studied by Berg et al., we evaluate the 

dimensionless frequency r2=27r ^y^^!^^ . =7rfro, where ro= — g^ ^"^ is the de- 
cay rate for concentration fluctuations. We notice that more than 75% of 
the data of Berg et al. are in the range Q<1. In the range where Q is small. 
The ratio R= J"'^}'^'^^ is linear in fi. This ratio is the most direct probe of 
the viscoelastic effect and one needs to concentrate on it. It is clear that for 

n«i 

3+Xn '0 

The result of plotting R vs. ^2 is shown in Fig.l. The one loop theory 
is shown by the dashed line. The slope of the data almost double. This is 
the discrepancy reported by Berg et al. In view of this differency, we have 
undertaken a two loop calculation of the frequency scale. In the low frequency 
end, we find a significant correction to the scale. This results in the solid 
line shown in Fig.l. In the high frequency end, the correction to the scale 
is similar. The data in the high frequency end is sparse. As far as we 



can tell, the scale that can be extracted from the data (a/ijis significantly 
smaller the low frequency result (aig). This is consistent with the calculation. 
The remaining difference between the experimental slope at low frequencies 
and the calculation can be attributed to the loops left out .The fact that 
including two loop correction is a significant effect, is a clear out pointer 
to the importance of higher order terms in perturbation theory at this end. 
We now outline the calculation involved. It should be noted that in the 
self consistent scheme that we employ, the two loop graphs corresponding 
to self energy insertions in the propagating lines have already been taken 
into account at the dressed one loop level. Consequently,treating the vertex 
correction alone enables us to provide a complete two loop order calculation. 

The two loop contribution to the shear viscosity r] (see Fig. 2) is given by 
the vertex correction diagram and can be written as 



772 fc^ 



1 f d^p f d^p [p^ - {k - p)^][q^ - {k - q)^] 



D-lJ {2n)D J (27r)^r7o(p2 + /t2)(g2 + ^2)|p'+^_^|2 

{PaTap{k)qp){paTap{k - p- q)qp) 

[-iuo + r(p, k) + r(A; - p, K)][-iuo + r(g, k) + r(A; - g, k) 



(10) 



In the above Tap^k) is the projection operator 6ai3 — ^r^and T{k,K) is the 

fully dressed order parameter relaxation rate. From the right hand side, we 

need to extract the 0(A;^2) term. We also need to average over all possible 

directions of k. 

Accordingly, 

p^-{k — p)"^ ~ -2k.p 

q^-{k- qf' ^ -2k.q 

and 

< {k.p){k.q)paTai3{k)qp >=^^^{SfJr^ 

Everywhere else in the right hand side, we may set k=0 in the right hand 
side of Eq.(lO). Thus the directional average oipaTapik — p — q)qp becomes 

< PaTai3{p + (l)lii > ^ud cau be written as 

T) Q SXTl 

< Parap{p + q)qa >= — ,^ ^2 (11) 

We can now write Eq.(lO) as 
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Specializing to D=3,we can replace the relaxation rate r(A;, k) by an accurate 
approximate to the full Kawasaki function as r{k,K) = rQk'^{k'^ + k,'^)'^ ,and 
we have 



d^p f d^q p'^q'^sin^9[l — 3cos^^] 



30r/ororoy (27r)3y {2T,f{p^ + l){q^ + l){p + q)^ 
1 



r]2 



[-1^ + pY(1 + p^)][-i^ + gY(l + g2)] 

The two loop contribution /nS'2 to InS is accordingly 

47r^ r d^p r d^q p'^q'^sin'^6[l — Scos'^O] 



(13) 



VoToJ {2nfJ {2Tif{l+p^){l + q'){p + q)^ 
1 
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[-1^ + pY(1 + p2)] [_^1] + g2^(l + ^2)] p2^(l + P^)q^^J{l + g2) 

(14) 

Our first task is to ensure that the zero frequency integrals reproduce the cor- 
rect two loop viscosity exponent Xn- Accordingly, from Equations(6)and(13), 
we find (the product tiqTq is fixed by the diffusion coefficient diagrammatics 
to be j^) 

Q 8 A 

V = r7i(/€, 0) + r]2iK, 0) ~ ^0T^^(1 + :r^)ln- + •••• (15) 

The coefficient of rjoln- in the above equation is within 2 % (the correction 
coming from two loop self energy insertion graphs and the dissipative four 
point coupling) is the exponent a;^. We can now carry out a Taylor expansion 
of the right hand side of Eq.(14). This yields the two loop contribution 02 to 
the scale factor aiQ as 

16 r d^p r d^qsin'^6{l — 3cos'^6)q'^ , . 

""'" isi (27r)3(l+p2)2 7 (27r)3(l + g2)|(p-'+5-)4 ^ ^ 

which yields 

«'o = 3[^ + A X ^ X 0.115] = 0.78 (17) 

lb OTT O 

Using the above scale, the solid line has been drawn in Fig.l. While this 
is a significant improvement over the one loop answer (dashed line), we still 
have an infinite number of loops left out and the combined effect should be 
to remove the remaining discrepancy. 

Turning to the high frequency side, evaluation of the appropriate integrals 
(Eq.l3) show that the behaviour of zero frequency viscosity to two loop order 

is 

8 , 8 ,„ A (I - /n2) + ^ X 0.85, 



while the high frequency hmit is 

8 , 8 A A^ A X 0.90, , , 

From Eqs.(18)and(19), the two loop a/i-is 

ah, = 0.23 (20) 

which is a 50% increase over the one loop result. 

Thus,we see that the scale factor associated with the frequency dependent 
viscosity near the critical point undergoes a significant enhancement both in 
the high and low frequency range, when the perturbation theory is carried 
out to two loop order. This helps us the reported discrepancy between the 
measurements of Berg et al.and the one loop calculation. 
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FIGURE - CAPTIONS 



Fig.l : The ratio of imaginary to real parts of the finite frequency shear 
viscosity plotted against the scaled frequency Q = -^^ for low values of 

Vt. The dashed line shows the one loop calculation, while the solid line 



is the two loop result. 



Fig. 2 : The diagrammatic expansion of shear viscosity to two loop order. 



